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The isotope elTects for Fe based superconductors are considered by including the phonon and 
magnetic fluctuations within the two band Eliashberg theory. We show that the recently observed 
inverse isotope effects of Fe, ape ~ — O.fS ± 0.03, [3| as well as the large positive isotope exponent 
(a « 0.35) can naturally arise for the magnetically induced sign revered s-wave pairing state within 
reasonable parameter range. Either experimental report can not be discarded from the present 
analysis based on the parameter values they require. The inverse and positive isotope effects mean, 
respectively, the interband and intraband dominant eletron-phonon interaction. We first make 
our points based on the analytic result from the square well potential model and present explicit 
numerical calculations of the two band Eliashberg theory. 

PACS numbers: PACS: 74.25. Kc, 74.20.Mn, 74.20.Rp 



Introduction - The alluring prospect of opening a key 
window to understanding the mechanism of high tem- 
perature superconductivity has attracted fierce research 
activities in the iron based pnictides.jl, 0, 0| The first 
step towards this path is to estabUsh the pairing symme- 
try and dynamics of the FeAs superconductors. Among 
the many ideas put forwardjparticularly appealing is the 
sign reversed pairing state. [1, @| It is the ground state of 
a two band superconductivity where both pairing order 
parameters, Ai and A2, on the two bands have full gaps 
while acquiring the tt phase shift between them. A re- 
pulsive interband interaction is turned to induce pairing 
by generating the sign reversal between the two order pa- 
rameters. The subsequent theoretical studies supported 
the idea of the sign reversed full gaps. The s± or stt 
pairing state seems to be able explain the experimen- 
tal observations indicating the full gap as well as a gap 
with nodes. 0, H, & [lO| The evidence for the full gaps 
came from the several independent ARPES experiments 
on single crystals [ll|, [l^, [ij] ■ The phase shift of tt may 
account for the absence of the Hebel-Slichter peak and 
the low temperature behavior of ~ in the spin-lattice 
relaxation rate 1/Ti,[l4, 15, 1^ which seemed difficult to 
understand in terms of a full gap scenario. 

As for the pairing interaction, one widespread school 
of thought maintains that magnetic fluctuations are inti- 
mately involved for the superconductivity. The supercon- 
ductivity emerges out of the interplay between the den- 
sity of charge carriers and strength of antiferromagnetic 
(AF) interaction such that their combined effect becomes 
optimum. This view seems natural from the overall phase 



diagram of the pnictides in the temperature and doping 
plane. The supercondcutivity emerges from the parent 
metallic antiferromagnetic state as the static AF order 
is suppressed and charge carriers are introduced through 
the chemical doping, [l'^ This view is further supported 
by the neutron scattering investigations showing a res- 
onance at a wave vector related to the AF order coin- 
cident with the onset of superconductivity. Within 
this school of thought combined with the theoretical es- 
timates of the small electron-phonon coupling constant, 
it came rather as a surprise that a large isotope coefficient 
a ss 0.34 — 0.37 was observed by Liu et al. with a sub- 
stitution of 5^Fe by ^''Fe in SmFeAsOi-^rF^, {x = 0.15) 
and Bai_:rKj;Fe2As2 (x — 0.4).[l9| Theoretical calcula- 
tion by Boeri et al. reported a very weak electron-phonon 
coupling constant {Xph « 0.21),[20j which is too small to 
account for the observed of the pnictides. Recall a is 
indeed very small for the cuprate supercondcutors par- 
ticularly near the optimal doping. [2l] | Perhaps even more 
surprising is the observation of the sizable inverse iso- 
tope effect of ape = —0.18 ± 0.03 by Shirage et al. in 
Bai_:,K^Fe2As2 (x = 0.4) fl] 

In view of this controversy, it is important that the Liu 
et al. and Shirage et al. results be confirmed. What we 
study here is to check if these experimental results are 
compatible with the current understanding of the pnic- 
tide superconductivity. We can estimate the parameters 
like the coupling constants that these experiments re- 
quire. These motivate us to revisit the issue of electron- 
phonon coupling in the Fe pnictides. We argue that the 
inverse isotope effect as well as the large positive effect 
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can arise with a reasonable parameter range within the 
widely held view of the magnetically mediated tt phase 
shifted s-wave pairing scheme. We first make our points 
based on the analytic result from the square well poten- 
tial model applied to the sign changing s-wave state, and 
substantiate them by the explicit numerical calculations 
of the two band Eliashberg theory. 

Idea and Model - As sketched below in the Formalism 
section, a straightforward application of the square well 
potential model to the case where there are AF spin fluc- 
tuations and phonons yields that the critical temperature 
Tc and isotope exponent a are given as follows: 
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where ujaf and ujph are the AF fluctuation and phonon 
frequencies, and A^j^^^/j is given by the sum/difference 
between the interband and intraband interactions as 



± (mier) , ^(intra) 

^AF/ph ~ ^AF/ph ^ ^AFjph ' 



(4) 



p/i ^ or < 0, the isotope effect can show 



Note that the Tc formula of Eq. ([T]) is reduced to the 
formula in Shirage et aLJJ in the limit that — 
and the quasi-particle renormalization is neglected. 
Also similar formulas were obtained before in different 
contexts. m ^ The result of Eqs. ^ and (HD, which 
will be backed up below by numerically solving the two 
band Eliashberg theory, explicitly shows that depending 
on whether A, 

the inverse or conventional behavior 

Two factors are behind this interesting effect: the — 
sign between the A^/f and A~^ of Eq. ([2]) and the pos- 
sibility of both signs of A~^ of Eq. The first factor 
comes from that the phonon and spin fluctuations enter 
the pairing kernel a|^ with the opposite sign as in Eqs. 
^ and (O, and the second factor from the sign change 
between the two gaps as suggested by the s± model. 
An inspection of a of Eq. ^ then reveals that the inverse 
isotope effect implies that the interband electron-phonon 
coupling is comparable with or stronger than the intra- 



band coupling. Note also that Ap'^*'^'' 



ph 



inverse isotope effect except when both are equal to zero. 
Also noteworthy is that when Ae// — A^^ « 0, Tc be- 
comes very small and a giant isotope effect will appear, 
|a| :$> 1/2. It seems that both the large isotope effect and 
the inverse isotope effect can arise in the magnetically in- 
duced TT phase shifted pairing for the pnictides within a 
reasonable parameter range. Either experimental report 



can not be discarded from the present analysis based on 
the parameter values they require. 

Although this discussion demonstrates the essential 
physics, it neglects the more general spin and phonon 
interactions and involves the square well potential model 
which is not a controled approximation. We will therefore 
present the results of the numerical calculations to show 
that inclusion of these effects still support this conclusion. 
We will show this by two band Eliashberg calculations as 
presented in what follows. 

Formalism - The two-band Eliashberg equation with 
an isotropic gap in each band may be written in the Mat- 
subara frequency as 



Ziiipn) = 1 
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where i,j = 1,2 are the band indices, aj„ = Pn — Pm, 
T the temperature, and p„ and pm are the Matsubara 
frequencies. The kernel A^^"* are given by 
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We take the phonon and spin fluctuations as 
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where F{e) is a truncated Lorentzian centered at uj^ 
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We took Tc = LOv, r — rc/2, and i? is a normalization 
constant such that deF{e) = 1. 24| The Tc — was 
chosen such that F{e) ^ to for small uj to simulate the 
overdamped AF spin fluctuations. 

Eqs. (O and ([S]) may be linearized near T ^ T^ and 
may be written as the following matrix form in the basis 

of Ai{ipn). 



Ai(ipn) = Mij{n, m)Aj{ipm), 

j m 

where the matrix M is defined as 
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The Tc is then determined to be the temperature at which 
the largest eigenvalue of the matrix M becomes 1. The 
isotope coefficient is given by 

1 d\nTc 



2 9 In ujpfi 
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Preliminary Analysis - Let us first consider some 
preliminary analytic analysis based on the square well 
potential model applied to the s± state. After neglecting 
the frequency dependence of x[^\iujn) and taking the 
standard approximation of Ai (w) = for < < ujph 
and Aib for Wph < OJ < ujaf, and similarly for Zi{uj), we 
can rewrite Eqs. ([TT|) and (fTU)) as 
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where the band indices take i,j — 1, 2, and 
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A straightforward manipulation taking advantage of the 
relations Aia = —A2a and Aib — —A2b of the s± 
state gives Eqs. ([T]), ([2]), and Q given above. Here, 

we used the definitions X^,^,^^ = xf^^^'^ ± xff^^'^ and 



A, 



AF/ph 
21 



and 



^AF/ph 

considered the simple case of xf^^^'^ 

^ AF/ph _ ^ AF/ph 
-^11 — -^22 

Before we present the results from the numerical com- 
putations of the linearized Eliashberg equation of Eq. 
(llOp . we first consider the approximate solutions from 
Eqs. (HI) and ^ to get ourselves oriented to the relevant 
parameter space. Let us first consider the inverse iso- 
tope effect in the context of the stt pairing state. This 
may be realized when the interband dominates the in- 
traband phonon interaction as discussed above. We take 
the simple parameterization of the interband interactions 



only, A^^/p^ - XAF/ph iKi 
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X 



AF/ph 
22 



0). From 



Tc « 3.0 meV and a ~ —0.2 reported by Shirage et al. for 
Bai_a;Ka;Fe2As2, we get Xaf = 1-53 and Xph = 0.189 for 
(jjph = i^AF = 20 meV, and Xaf = 0.944 and Xph = 0.146 
for Wph — 20, uaf — 30 meV. The obtained Xaf > Xph 
is consistent with the magnetically induced pairing. 

We now turn to the a > case, which one of the au- 
thors considered previously. [2^ This may be realized in 
the context of the stt pairing state by X^p = Xaf and 
A^^ = ±Xph- From Tc ~ 3.5 meV and a « 0.35 reported 
by Liu et al. for SmFeAsOi_a;Fa;, we get Xaf ~ 0.836, 
and Xph — 0.738 by solving Eqs. H]) and ([2]) for iVph = 20 
meV and loaf — 10 meV. The solution for uiph — uiaf — 
20 meV does not satisfy Xaf > Xph relevant for the mag- 
netically mediated superconductivity. The inverse or pos- 
itive isotope effect is naturally obtained when luaf > ^ph 



or <. This was pointed out by Shimahara in a different 
context. 122] 

In the derivation of Eqs. H]) and ([2]), we assumed 
UJAF > ^ph', if UJAF < i^ph, they are not valid. Also 
they are based on the uncontrolled double well approx- 
imation. The correctness of the results is confirmed by 
numerical calculations we are presenting below. 
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FIG. 1: The transition temperature Tc and isotope exponent 
Of as a function of the phonon frequency ujph- ujaf ~ 30 meV, 
Xaf ~ 2.26, Xph = 0.5, and all the intraband couplings are 
set equal to 0. The solid lines are the plots of Eqs. ([iJ and 
([2|). See the text for their parameters. 

Numerical Results - In Fig. 1, we show the results of 
numerical calculations of the inverse isotope effect. We 
took, guided by the above discussion, ujaf = 30 meV, 
and adjusted the interband Xaf and Xph with all intra- 
band interactions set to 0, so that the and a give 
the right values for Bai_a;Ka;Fe2As2. With A^^ = 2.26 
and Xph = 0.5, ujph was varied between 10 and 50 meV. 
The solid squares and triangles show the a and Tc, re- 
spectively. The lines are the plots from Eqs. ^ and 
(121). Xaf = 0.876 and Xph = 0.104 were selected such 
that Tc and a agree with the numerical calculations at 
ujph ~ 20 and ujaf = 30 meV. The interaction param- 
eters are somewhat different between the numerical and 
(approximate) analytic calculations. 

We then consider the positive isotope effect and present 
the results in Fig. 2. We took ujaf = 15 meV, interband 



A 



AF 



1.06, intraband A 



ph 



0.950, and varied ujph be- 
tween 10 and 50 meV as before, with other interactions 
set to 0. The solid triangles and squares show the a and 
Tc, respectively, and the lines are from the analytic ex- 
pression with Xaf = 0.878 and Xph = 0.491 to agree with 
the numerical results at ujaf = 15 and ujph — 20 meV. 

Although we have not fine-tuned the parameters 
to match the observed Tc and a exactly for the 
Bai_2:K:rFe2As2 and SmFeAs0i_3;Fa;, some features are 
worth noting. First, the parameters for the numerical 
and analytical calculations to give the observed Tc and 
a were somewhat different. This is because Xaf and 
Xph vary sensitively as ujaf is varied (but not so sensi- 
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FIG. 2: Tc and a as a function of uiph- i^AF ~ 15 meV, 
\af = 1.06, A^^"'™' = 0.950. The solid lines are the plots of 
Eqs. (O and ([2}. See the text for their parameters. 



lively on LUph) in this region of parameter space as one 
can see easily by solving Eqs. ^ and ([2]). Second, it 
is necessary that the phonons are predominantly of in- 
terband/intraband character for the inverse/positive iso- 
tope effects. And, third, the phonon coupling constant 
somewhat larger than the theoretical calculation seems 
necessarv. [20I] But, this might be understood by including 
other effects like the different density of states between 
the bands. [2^ 

Remarks - We have shown that both the inverse iso- 
tope effect and large isotope exponent reported in the 
pnictides may arise in the magnetically induced tt phase 
shifted pairing with reasonable parameters. Either ex- 
perimental report can not be discarded from the present 
analysis based on the parameter values they require. The 
dominant interband (intraband) phonon interaction pro- 
duces the inverse (positive) isotope effect. Note that the 
contradicting inverse and positive isotope effects were ob- 
served in the same compound, Bai_a;Ka;Fe2As2 {x — 0.4), 
by Fe substitution. The Liu result was obtained on the 
sample from the solid state reaction whereas the Shirage 
result from the high pressure technique. If either result 
turns out to be reproducible, it will point to the interband 
or intraband dominant electron-phonon interaction in the 
pnictides. Liu et al also observed that O isotope effect 
is quite small, ao ~ —0.06. In the framework proposed 
here, this may be understood in terms of A^^ « « 
for the O ions. It remains to be checked if this antici- 
pation is correct. It will be also interesting to see if the 
spin density wave isotope exponent reported by Liu et al 
may be understood within the proposed idea here. 
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